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Abstract 

The dispersion relations of carbon nanotubes are often obtained cross-sectioning those of graphene 
(zone-folding technique) in a rectangular region of the reciprocal space, where it is easier to fold 
the resulting relations into the nanotube Brillouin zone. We propose a particular choice of the 
unit vectors for the graphene lattice, which consists of the symmetry vector and the translational 
vector of the considered carbon nanotube. Due to the properties of the corresponding unit vectors 
in the reciprocal space, this choice is particularly useful for understanding the relationship between 
the rectangular region where the folding procedure is most easily applied and the overall graphene 
reciprocal space. Such a choice allows one to find, from any graphene wave vector, the equivalent 
one inside the rectangular region in a computationally inexpensive way. As an example, we show 
how the use of these unit vectors makes it easy to limit the computation to the bands nearest to 
the energy maxima and minima when determining the nanotube dispersion relations from those of 
graphene with the zone-folding technique. 
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1. INTRODUCTION 



Carbon nanotubes are cylindrical structures with diameters that are usually in the few 
nanometer range and lengths up to tens of microns. Due to their high mechanical strength 
and thermal conductivity and to their unusual electronic properties, carbon nanotubes con- 
stitute a very promising material for many applications [2j , such as active devices, intra- 
chip interconnections, field emitters, antennas, sensors, scanning probes, reinforcement for 
composite materials, energy and hydrogen storage. In particular, from the electronic point 
of view, they can behave as metallic or semiconducting materials, depending on their geo- 
metrical properties ^, 0, Q] . 

A single-wall carbon nanotube can be described as a graphene sheet rolled, along one of 
its lattice vectors (the so-called chiral vector), into a cylindrical shape. As a consequence 
of the closure boundary condition along the chiral vector, only a subset of graphene wave 
vectors, located on parallel lines, are allowed. Therefore the dispersion relations of carbon 
nanotubes are often found cross-sectioning those of graphene along such lines (zone-folding 
technique) . The cross- sect ions are usually taken in a particular rectangular region of 
the graphene reciprocal space (which can be seen as a primitive unit cell of the graphene 
reciprocal lattice). 

Here we introduce an unusual choice of the unit vectors in the graphene direct and recipro- 
cal space, which, as a result of a direct geometrical relation with such a rectangle, makes it 
clearer how the overall reciprocal space can be obtained by replicating the rectangle. 
This allows a more complete understanding of the results of the zone-folding method, ex- 
plaining, for example, how the periodicity of the nanotube energy bands arises from the 
computational procedure. 

We also show how such unit vectors make it easy to find, from any point of interest of the 
graphene reciprocal space, the equivalent wave vector located inside the above-mentioned 
rectangular region. At the end of the present communication, we apply this procedure to 
the graphene degeneration points, in order to compute just the nanotube bands nearest to 
the energy maxima and minima. 

In Figs. [H and [2] we show the graphene lattice in the direct and reciprocal space, respec- 
tively, and the (x, y) reference frame that we have used in the following. The graphene 
lattice structure in the real space can be seen as the replication of the graphene rhomboidal 
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unit cell shown in Fig. [T] (containing two inequivalent carbon atoms A and B) through linear 
combinations with integer coefficients of the lattice unit vectors a\ = (\/ r 3a/2)x + (a/2)y 
and a 2 = (V3a/2)x — (a/2)y. Correspondingly, the unit vectors of the graphene reciprocal 
lattice are b\ = (2n / (^/3a))x + (2ii/a)y and b 2 = (27i/(y/3a))x — (2ii/a)y, which are reported 
in Fig. El along with the graphene hexagonal Brillouin zone. 

An (n, m) carbon nanotube is obtained rolling up a graphene sheet along its chiral vector 
Ch = nSi + ma 2 ; the circumference of the nanotube is consequently equal to the length of 
this vector: L = \Ch\ = a\/ n 2 + m 2 + nm. 

If we define dpt as the greatest common divisor of 2m + n and 2n + m, we have that the 
lattice unit vector of the nanotube (which represents a ID lattice) is the so-called trans- 
lational vector T = t\di + t 2 a 2 of the unrolled graphene sheet, parallel to the nanotube 
axis and orthogonal to C^, where t\ and t 2 are relatively prime integer numbers given by 
ti = (2m + n)/dn and t 2 = — (2n + m)/d R . Therefore, the rectangle having as edges the 
chiral vector Ch and the translational vector T represents the unit cell of the nanotube, 
which repeats identically along the nanotube axis with a lattice unit vector T, the length 
of which is equal to T = \T\ = y/3L/dn. The number of graphene unit cells inside the 
nanotube unit cell is equal to iV = 2L 2 /(a 2 dn). 

The coordinates of all the N points identifying the graphene unit cells inside the rectangular 
region representing the nanotube unit cell in the unrolled graphene sheet are defined (apart 
from translations by an integer number of Ch and T vectors) by integer multiples of the 
so-called symmetry vector R = pa\ + qa.2, where p and q are two relatively prime integer 
numbers, univocally determined by the two relations: t\q — t 2 p = 1 and < M < N (where 
we define the quantity M = mp — nq). In particular, we have that NR = Ch + MT [6]. 
In Fig. [T]we show all of these vectors in the direct space for the nanotube (10, 0) (for which 
C h = 10ai, L = 10a, d R = 10, f = a x - 2a 2 , T = y/3a, N = 20, R = a x - a 2 and M = 10). 
Since in the direct space the nanotube is a one-dimensional lattice with a unit vector T and 
with a T wide unit cell along the nanotube axis, in the reciprocal space its unit vector is equal 
to K 2 = (2n/T)T = (mbi — nb 2 )/N and its Brillouin zone is represented by the values of 
the wave vector k (along the nanotube axis) which satisfy the inequality —n/T < k < 7r/T. 
Many physical properties of the nanotubes, such as the energy dispersion relations, can be 
found from the corresponding quantities of graphene using the zone-folding technique 
Indeed, as a consequence of the closure of the graphene sheet to form the carbon nanotube, 
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FIG. 1: The described quantities in the direct space for a (10,0) nanotube. The heavy-shaded 
rhombus is the graphene unit cell, while the light-shaded rectangle represents (once the graphene 
sheet has been rolled up) the nanotube unit cell. 

we have to enforce that the graphene electron wave function e tk ' r u(k) (where u(k) is a Bloch 
lattice function) has identical values in any pair of points r and r + Ch- To satisfy the 
resulting relation e lk ' Ch = 1, the component along C h of the wave vector has to be equal to 
an integer multiple of the vector K\ = (2ir/L)C h = (— tS\ +tib 2 )/N. If we cross-section the 
graphene dispersion relations in correspondence of the parallel lines (separated by a distance 
\Ki\ = 27r/L) containing the allowed graphene wave vectors and we fold such sections into 
the nanotube Brillouin zone, we find the relations for the carbon nanotube. This procedure is 
applied to a region of the graphene reciprocal space containing all and only the inequivalent 

— * 

graphene wave vectors. In particular, the rectangle having as edges the vectors NKi and 

— # 

K 2 has these properties (as we have demonstrated in the Supplementary Information) and 
corresponds to the region that is usually implicitly chosen to apply the zone-folding method, 
because here, considering only the allowed graphene wave vectors, we obtain N segments 
with a width equal to 2ir/T, that can be folded into the Brillouin zone of the nanotube 
by simply taking the component along T of each graphene wave vector (such a component 
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FIG. 2: The described quantities in the reciprocal space for a (10,0) nanotube. The light-shaded 
hexagon is the graphene Brillouin zone, while the heavy-shaded rectangle is the rectangular area 
specified by the vectors NK\ and K 2 . 

becomes the nanotube wave vector). 

In Fig. [2] we show the quantities in the reciprocal space for the nanotube (10, 0) (for which 
K x = (27r/(10a))Cft = (26i + b 2 )/20 and K 2 = (27r/( v / 3a))f = -b 2 /2). 
We propose an alternative choice of the graphene unit vectors that allows to clarify the 
relation between the previously described rectangular region specified by the vectors NK\ 
and K 2 and the overall graphene reciprocal space; in particular it makes it easy to find, for 
any given wave vector, the equivalent wave vector belonging to such a rectangular area. As 
we shall see, this can be very useful when we apply the zone-folding method cross-sectioning 
this rectangle. 
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2. ALTERNATIVE CHOICE OF THE GRAPHENE UNIT VECTORS 



Let us consider the graphene sheet forming, once rolled up, a carbon nanotube with chiral 
vector Ch, translational vector T and symmetry vector R. We propose, as an alternative 

— # 

choice of the graphene unit vectors in the direct space, exactly the pair of vectors cIa = R 
and Sb = T: 

a A = R = pa 1 + qa 2 = ^Vs(p + q)x + - (p - q)y (1) 

a B = f = tiai + t 2 a 2 = |\/3(ti + h)x + |(*i - t 2 )y • (2) 

To verify that these two vectors can actually be used as graphene unit vectors in the direct 
space, we have to demonstrate that their linear combinations with integer coefficients yield all 
and only the lattice vectors of such a space, i.e. the vectors that are also linear combinations 
with integer coefficients of at and a 2 . 

Since both R and T are linear combinations with integer coefficients of Si and a 2 , every 
linear combination with integer coefficients of R and T is also a linear combination with 
integer coefficients of a*i and a 2 . 

On the other hand, in order to demonstrate that every linear combination with integer 
coefficients of a\ and a 2 is also a linear combination with integer coefficients of R and T, it 
is useful to consider the linear system consisting of the two following known relations: 

pax + qa 2 = R 
Mi + t 2 a 2 = T 
Solving for a\ and a 2 , we find that 

(-t 2 p + qti)3i = -t 2 R + qf . \a 1 = -t 2 R + qf 

or equivalently < (4) 

(t^ - pt 2 )a 2 = t 1 R - pf I a 2 = t 1 R -pT 

(recalling that tig — — 1)? which means that cL\ and 02 are linear combinations with 
integer coefficients of R and T. Consequently, also every linear combination with integer 
coefficients of a% and a 2 is a linear combination with integer coefficients of R and T. 
Using the coordinates of R and T in the (x, y) reference frame of Fig. (TJ and introducing a 
unit vector z that is orthogonal to the plane x, y of the graphene sheet and forms a right-hand 
reference frame (x, y, z) with x and y, we have that 

fxz = ~(t 2 -t 1 )x-^VS(t 2 +t 1 )y (5) 



z X R = ~(q-p)x + -Vs(q + p)y (6) 



i2.(Tx«) = -^-(* 1 g-t a p) = -^- , (7) 
where we have used again the relation t\q — t 2 p = 1. Therefore, the corresponding unit 



vectors of graphene in the reciprocal space are (using t 



unit vectors in the direct and in the reciprocal space js]): 

r a B x z T x z 2n , . 2 it . 

b A = 2vr- — ^ = 2tt- - = -= (t 2 - h)x (t 2 + h)y 

a A -(a B xz) R-(Tx z) aV3 a 

( 2tt a 2tt \ / 2tt a 2tt \ -> r 

= ~h —t=-x H y + ti -7=^z y = -fe&i + tw = 

Vv3a a / \V3a a / 



re well-known relations between the 



-t 2 (nKi +tiK 2 )+ tiimKt + t 2 K 2 ) = (-t 2 n + t 1 m)K 



NKt = N^-C h (8) 



and 



-> z x a a z x R 27r 2n 

b B = 2n- — — = 2n— - = — =(o — p)x H (q + p)y = 

a A -{a B xz) R-(Tx z) ay/T 1 V) a W F,y 

( 2vr 2tt \ / 2vr 2tt \ r r 

= q{nKi + t x K 2 ) - p{mK x + t 2 K 2 ) = (qn - pm)K 1 + {qt x - pt 2 )K 2 = 

2tx - 2n ~ 

= -MK l +K 2 = -M—C h + —T , (9) 

— * —* 

which are linear combinations with integer coefficients of the vectors K 1 and K 2 and have 
components along C h and T: b Ac = N{2n/ L) , b Ar = , b Bc = —M(2ii/L) , b Br = 2n /T . 
To find the previous results we have used, besides the fact that tiq—t 2 p = 1 and M = mp—nq, 
the relations: — t 2 n + t\m = N, b\ = nK\ + t\K 2 and b 2 = mK\ + t 2 K 2 . 
Indeed, using the relations listed in the Introduction, we obtain that 

2m + n 2n + m 2 . 2 2 . 2 / L\ 2 . , 

mbi — nt 2 = m — h n — = — (m + nm + n ) = — I — I — N . (10) 
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The relations between the vectors b\ and b 2 and the vectors K 1 and K± can instead be found 
starting from the identities: 



K 2 = j-fimbi ~ n h) 
and solving for bi and b 2 this system of two equations. We find that 

iratx - nt 2 )b x = N{nK\ + hK 2 ) 
(mti - nt 2 )b 2 = NimKx + t 2 K 2 ) 

and thus (using the fact that mt\ — nt 2 = N) that 

2tt - 2vr - 

61 = nKi + tiK 2 = n—C h + h—T 

2tt * 2tt * 

6 2 = mifi + t 2 K 2 = m—C h + t 2 —T 



(11) 



(12) 



(13) 



From this result, it is also apparent that the components of bi and b 2 along Ch and t are: 
6i c = n(27r/L) , &i T = ti(2vr/T) , 6 2c = m(27r/L) , 6 2t = t 2 (2n/T) . 

3. APPLICATIONS 

This choice of unit vectors helps us to understand the relation between the rectangular region 
having as edges the vectors NKi and K 2 and the overall graphene reciprocal space. Indeed, 
such a rectangular region contains all and only the inequivalent graphene wave vectors and 
can therefore be considered as a primitive unit cell of the graphene reciprocal lattice. In 
particular, considering as unit vectors of the graphene reciprocal space &a and &b (which 
have a clear geometrical relation with the considered region), we have that the overall re- 
ciprocal space can be spanned translating the rectangular region by vectors that are linear 
combinations, with integer coefficients, of &a (which is exactly equal to NKi, the base of the 
rectangle) and 6# (which has a component along T equal to \K 2 \, the height of the rectangle, 
and a component along Ch equal to —M\Ki\, i.e. an integer number of times the distance 
\K\\ between the segments along which we take the N cross-sections inside the rectangle in 
the zone-folding method). Therefore, the overall graphene reciprocal space is spanned by 
rows (parallel to Ch) of equivalent rectangles, with each row shifted along Ch by —MK\ 
with respect to the adjacent one (as we show in gray in Fig. [2] for the particular case of a 
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(10,0) nanotube, where M = 10, N = 20 and therefore N = 2M). 

This clarifies the result obtained by cross-sectioning the graphene dispersion relations along 
the N parallel lines (separated by a distance |i?i|) to which the N parallel segments used in 
the zone-folding method belong. Since the parallel rows of rectangles spanning the graphene 
reciprocal space have a nonzero relative shift along Ch (and therefore the generic graphene 
wave vectors k and k + K 2 are not equivalent), the relation obtained from each single cross- 
section in general is not periodic with period equal to \K 2 \ (the width of the nanotube 
Brillouin zone). Nevertheless, since the relative shift along C h between rows of rectangles is 
an integer multiple of the distance \K\\ between the N parallel lines, we find, starting from 
a wave vector k on one of the N lines and moving by \K 2 \ along T and by a proper multiple 
of \Ki\ along Ch, a wave vector equivalent to k on another of the N lines of allowed wave 
vectors. Therefore the overall set of relations, obtained drawing all the N cross-sections on 

— * 

the same one-dimensional domain, is indeed periodic with period \K 2 \. This is in agreement 
with the fact that the resulting relations are the nanotube dispersion relations, that have to 
be periodic with a period equal to the width of the nanotube Brillouin zone. 
The proposed alternative choice of graphene unit vectors is particularly useful for the de- 
termination, for any given wave vector, of the equivalent wave vector within the previously 
discussed rectangular region of the graphene reciprocal space. 

— * — * — * 

Indeed, given a graphene wave vector k, if we use b A and b B as unit vectors in the reciprocal 
space, all the wave vectors equivalent to k can be written as k eq = k + ab A + fibs, with a 
and (3 integer numbers. Thus the corresponding components k^ and along Ch and T, 
respectively, are: 

k e J = k c + ab Ac + (3b Bc = k c + aN— - (3M— 

Li Li 

14 

k^ = k T + ab Ar + l3b Br = k T + f3— 

— * 

Since we want to find the particular k eg belonging to the rectangular region, such components 
have to satisfy the following relations: 



< k e J < N— 



.2tt 

T 



T T ~ T 



or equivalently 



27T 

< k e S < N— 

. (15) 

0<k? + %<% 
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Substituting the expressions of and Wp into these inequalities, we find: 



2tx 



.2tt 



,2tt 



< [he- P M j;J + aN j; < N ~ 

7T\ 2vr 2tt 
< [k T + -) + (3— < — 



T 



(16) 



and thus 



.2tt 



L 



T ~ T 



2ty 



a N— ) <k c - (3M— < (-a + 1) TV— 



.2tt 



/2tt 



or equivalently 



)<fc T + ^<(-/3 + l) 
k c - /3M(2tt/L) 



2tt\ 

t 7 ; 



(17) 



-a < 



N(2tt/L) 



< -a + 1 



M 2tt/T ~ M 



(18) 



The values of a and /3, and consequently the vector k eq , can be easily found using the fact 
that the second inequality contains only f3. Indeed, from the second inequality we find that 
(using the ceiling and floor functions): 

~ k T + (tt/T) 



2n/T 

k T +{<K/T)\ 
2vr/T 



7T 



+ 1 if k T + — > 



(19) 



7T 



if k T + - < 



Once the value of (3 is found, the quantity kc — (3M(2n/L) in the first inequality is known 
and thus from first inequality we obtain that: 



a 



k c - PM(2ir/L) 
N(2tt/L) 

k c -pM(2ir/L)\ ' 



if k c - pM^- > 



if k 



c 



PM^- < 



(20) 



N(2tt/L) 

Clearly the second inequality of the systems (fTB"]) - (fl8"]) does not contain a only because, with 
our particular choice of unit vectors, %a has a zero component along T (i.e. bA T = 0). 
In the following we show an application of this procedure for the optimization of the zone- 
folding computation of the nanotube energy bands. 

As we have described, the nanotube dispersion relations can be obtained cross-sectioning the 
bands of graphene (computed for example with the tight-binding method), in the rectangular 
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region of the reciprocal space specified by the vectors NKi and K 2 , along the N equidistant 
segments, parallel to T, containing all the wave vectors of the region with component along 
Ch multiple of \K\\ = 2n/L. In this way, cross-sectioning the two energy bands (bonding 
and anti-bonding) of graphene, we obtain 2N dispersion relations that, once folded into the 
nanotube Brillouin zone (which coincides with the first segment, along the nanotube axis), 
form the nanotube energy bands. 

Among these 2N energy bands, the most interesting ones are the lowest conduction bands 
and the highest valence bands, because these are the regions where the charge carriers 
localize. These bands are obtained cross-sectioning the graphene dispersion relations near 
the particular graphene wave vectors 

l -(2b lc +b 2c )C h + l -{2b lT +b 2T )f , (21) 

~(6 lc + 2b 2c )C h + i(6a T + 2b 2T )f (22) 

and their equivalent wave vectors, where the graphene energy bands have their maxima 
and minima (and in particular are degenerate). Therefore we need to find the wave vectors 
equivalent to K and K' inside the rectangular region where we take the cross-sections of 
the graphene dispersion relations (we shall find just one wave vector equivalent to K and 
just one equivalent to K'). This is done applying the previously described procedure to K 
and K' whose components along Ch and T are given by Eqs. (I21|) - (j22j) . In particular, we 
know |6| that the components along T of the graphene wave vectors equivalent to K and 
K' and belonging to the rectangle can only assume the values or ±27r/(3T) and therefore 
are well inside the considered region, away from the boundaries. Once we have found these 
wave vectors, we can cross-section the graphene dispersion relations just along the segments 
in their proximity, if we want only the mentioned most relevant energy bands. This leads 
to a nonnegligible reduction of the computational effort. 

For example, in Fig. [3] we show the results obtained for a (10, 0) carbon nanotube (for 
which, as we have seen, N = 20). The graphene dispersion relations have been computed 
with the tight-binding method, in which we have considered only the 2p z orbital for each 
atom and have included the effect on each atom of up to the third-nearest neighbors. In this 
computation we have used, for the tight-binding parameters, the values found in [q| fitting 
in the optical range the results of an ab initio calculation performed with the SIESTA code: 



jt» 2tt a 2tt^ 1 r -.. 

K = — X +—y=~(2b 1 + b 2 ) 

V3a 3a 3 



2tt 2tt 1 r -♦, 
= V3a X ~3a V= 3 + 
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e 2p = -2.03 eV, 70 = -2.79 eV, s = 0.30, 71 = -0.68 eV, si = 0.046, 72 = -0.30 eV and 
S2 = 0.039 (with the notation used in 9j). We have considered 10001 wave vector values on 
each of the N parallel segments along which the graphene energy bands are cross-sectioned 
(in order to keep our code as general as possible, we did not exploit the particular simmetry 
properties of the achiral (10,0) nanotube). The curves shown in the figure represent all 
the 2N = 40 (partially degenerate) bands (20 bonding bands and 20 anti-bonding bands) 
of the (10, 0) nanotube. In particular, with the thin black lines we represent the 8 bands 
(degenerate in pairs) obtained cross-sectioning the graphene dispersion relations along the 
segments closest to the two graphene degeneration points inside the considered rectangle, 
and with the thick black lines we report the portions of these bands obtained taking the 
cross-sections only in the circular regions (with radius equal to (5/6) |i?i|) centered around 
the two graphene degeneration points. The data represented with the (thin and thick) black 
lines have been obtained with the previously described method. 

In our simulations we have found that, on a Pentium 4 at 2.4 GHz, the time spent to find all 
the bands of the (10, 0) nanotube is 200 ms, five times greater than the time (40 ms) taken 
by the modified version of the program, which computes only the bands closest to the two 
graphene degeneration points: the computational time substantially scales proportionally 
to the number of computed bands. Computing only the parts of such bands closest to the 
graphene degeneration points we have a further speed-up: in this case the time spent becomes 
10 ms. The proposed improvement becomes actually useful in the situations in which many 
calculations of this type need to be performed, leading to significant computational times. 
Incidentally, we note that an alternative method to calculate only the most relevant nanotube 
energy bands could consist in taking the cross-sections of the graphene dispersion relations 
(along the parallel lines corresponding to the allowed wave vectors) in the hexagonal Brillouin 
zone of graphene (which evidently contains all and only the inequivalent graphene wave 
vectors), instead of inside the considered rectangular region. In this case the positions 
of the maximum and minimum points are well known (K and K' are at the vertices of 
the hexagon) and thus the regions of interest are clearly located. Following this method, 
however, in order to fold the computed cross-sections into the nanotube Brillouin zone (the 
segment of the T axis characterized by — tt /T < k T < n/T), it is not sufficient to consider 
the kx component along T of the graphene wave vector (the absolute value of this component 
can be greater than 7r/T), but we also need to find the nanotube wave vector equivalent to 
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FIG. 3: Results of the proposed optimization of the tight-binding method in the case of a (10, 0) 
nanotube. The thin black lines are the (doubly degenerate) bands obtained cross-sectioning the 
graphene dispersion relations along the segments closest to the two graphene degeneration points 
inside the considered rectangle and the thick black lines are the parts of these bands obtained 
taking the cross-sections only in the circular regions (with radius equal to (5/6) | ) centered 
around the two graphene degeneration points. The gray curves represent the nanotube dispersion 
relations that with our procedure we do not need to compute, if we are not interested in them. 

hr inside the nanotube Brillouin zone. Moreover, since the graphene degeneration points 
are on the boundary of the sectioned hexagonal region, in order to obtain the nanotube 
bands around their minima and maxima we have to properly join the results computed 
cross-sectioning the graphene dispersion relations near the six vertices (each of which gives 
only parts of the desired bands); this strongly increases the algorithmic complexity of the 
involved computations. 

4. CONCLUSION 

We have proposed an alternative choice of unit vectors for the graphene sheet which, once 
rolled up, forms a carbon nanotube. These vectors, which depend on the considered nan- 
otube, are closely related to the rectangular region of the graphene reciprocal space where 
the zone-folding method is most easily applied, and allow us to better understand the rela- 
tion between this rectangular area and the overall reciprocal space. In particular, we have 
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shown that our choice of unit vectors can be exploited to find, from any graphene wave 
vector, the equivalent wave vector inside the rectangular region and can therefore be useful 
whenever the zone-folding technique is applied to obtain the physical properties of the car- 
bon nanotube from those of graphene. As an example, we have presented an application to 
the optimization of the tight-binding calculation of the carbon nanotube energy dispersion 
relations, with a significant reduction of computational times. 
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Supplementary Information 

In the following, we provide a demonstration of the fact that the rectangular region of the 

— * — * 

graphene reciprocal space defined by the vectors NK\ and K 2 contains all and only the 
inequivalent graphene wave vectors and can therefore be considered as a primitive unit cell 
of the graphene reciprocal lattice. 
This rectangular region has an area 

2vr 2tt _ 2L 2 2vr d R 8tt 2 

l ' Y~^d~ R L 7E~V^ ' (23J 

that is equal to \bi x b 2 \ and thus to the area of the graphene Brillouin zone. Thus our 
assertion is automatically verified if we demonstrate that such a region does not contain two 
distinct but equivalent (i.e. differing only for a linear combination with integer coefficients 
of b\ and b 2 ) graphene wave vectors k = k c Ch + k T f and k' = k + N 1 b 1 + N 2 b 2 = k' c C h + k' T T 
(with Ni and N 2 two integer numbers). 

— * 

Indeed, if k is inside the rectangular region, we have that 

< k c < N^- { —N^j- <-k c <0 

and therefore < " . (24) 

7T 7T i 7T 7T 

--<k T <- y--<-k T< - 
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If also k! is inside the rectangular region, we have that 



< k' r < N^- 



0<k c + NJtc + N 2 b 2c < N^- 



7T , . 7T 

-— < k' T < — 

T T — rp 



or equivalently 



7T 



T <k T + N 1 b lT + N 2 b 2T < - 



From these last relations, exploiting the inequalities for — k c and — k T , we have that 

27T 27T 27T 

—N— < -k c < N x b lc + N 2 b 2c < -k c + N— < N— 

Lj Lj Lj 

< -kr ~ t^< Nih T + N 2 b 2r < -k T + — < 



T 



This implies that that 



2vr 2tt 
—N— < + N 2 b 2c < N— 

Lj Lj 

27r lr , A7I 27T 

-— < iVx6 lT + N 2 b 2T < — 



T T 



and thus 



\Ntb lc + N 2 b 



2c 



2tt 



|JVi6i t +iV 2 6 2T | < — 



(25) 



(26) 



(27) 



(28) 



Substituting into (j27| the values of the components of b\ and 62 along Ch and T, we have 
that 

|JVin + N 2 m\ < N 
1^1 + N 2 t 2 1 < 1 

Being iViti + A^2^2 an integer number, the second inequality is equivalent to the relation 
Niti + iV 2 t2 = 0. Since N% and N 2 have to be integer numbers, and t\ and t 2 are two relative 
prime integer numbers, this identity is satisfied only if N\ = £t 2 and N 2 = —it\, with I an 
integer number. With these values of N\ and A^2, we have that Nin + N 2 m = i{t 2 n — t\vn) = 
—£N (as we have seen, mt\ — nt 2 = N). Thus the first inequality of (128]) becomes \£\N < N 
and is satisfied only if I = 0. This means that Ni = £t 2 = and N 2 = —£t± = and thus k' 
is identical to k, as we wanted to prove. 
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